Scalable quantum computation via local control of only two qubits 



0^ 
O 
O 
(N 



Oh 



> 
cn 
i> 
m 
rn 

o 

O 



X 




Daniel Burgarth^-^, Koji Maruyama^, Michael Murphy'^, Simone 

Montangero'^, Tommaso Calarco'^ Franco Nori^'^, and Martin B. Plenio^ 
^IMS and QOLS, Imperial College, London SW7 2PG, UK 
^Advanced Science Institute, The Institute of Physical and Chemical Research (RIKEN), Wako-shi, Saitama 351-0198, Japan 
^Institiit fr Quanteninformationsverarbeitung, Universitt Ulm, D-89069 Ulm, Germany 
'^ECT*; 38050 Villazzano (TN), Italy and 
^Physics Department, University of Michigan, Ann Arbor, Michigan, 48109, USA 

We apply quantum control techniques to control a large spin chain by only acting on two qubits at one of its 
ends, thereby implementing universal quantum computation by a combination of quantum gates on the latter 
and swap operations across the chain. It is shown that the control sequences can be computed and implemented 
efficiently. We discuss the application of these ideas to physical systems such as superconducting qubits in 
which full control of long chains is challenging. 
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Controlling quantum systems at will has been an aspira- 
tion for physicists for a long time. Achieving quantum con- 
trol not only clears the path towards a thorough understanding 
of quantum mechanics, but it also allows the exploration of 
novel devices whose functions are based on exotic quantum 
mechanical effects. Among others, the future success of quan- 
tum information processing depends largely on our ability to 
tame many-body quantum systems that are highly fragile. Al- 
though the progress of technology allows us to manipulate a 
small number of quanta quite well, controlling larger systems 
still represents a considerable challenge. Unless we over- 
come difficulties towards the control over large many-body 
systems, the benefits we can enjoy with the 'quantumness' 
will be severely limited. 

One encouraging result is that almost any coupled quantum 
system is controllable in principle, even by steering a single 
particle only Jj]]. The question of controllability in this type 
of situation can be described by the theory of quantum con- 
trol which uses Lie algebraic arguments. This is inter- 
esting from the theoretical point of view; however, can it be 
practically useful from the quantum computing perspective? 
Problems we need to contemplate before attempting to build 
a large quantum computer using quantum control are as fol- 
lows. First, the control criterion is generally not computable 
for large systems. Second, even if the question of controlla- 
bility can be answered positively for specific systems |3], the 
precise sequence of actual controls (or 'control pulses') are 
generally not computable. And third, even if they can be com- 
puted, the theory of control tells us nothing about the overall 
duration of the control pulses to achieve a given task, and it 
might take far too long to be practically relevant. 

The usual approach to circumvent these problems focuses 
on systems that are sufficiently small, so that we do not al- 
ready require a quantum computer to check their controllabil- 
ity and to design control pulses. In such a case, the theory of 
time optimal control [4] can be used to achieve impressive im- 
provements in terms of total time or type of pulses required in 
comparison with the standard gate model. More complicated 
desired operations on larger systems are then decomposed 
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Figure 1 : (color online) Our approach for universal quantum compu- 
tation works on a chain of spins. By modulating the magnetic field 
B\{t) on qubit 1, we induce information transfer and swap gates on 
the chain (red and green lines). The states of the qubits from the un- 
controlled register can be brought to the controlled part. There, the 
gates from a quantum algorithm are performed by local operations. 
Afterward, the (modified) states are swapped back into their original 
position. 



('compiled') into sequences of smaller ones. Yet, the feasibil- 
ity of this approach is ultimately limited by the power of our 
classical computers, therefore constrained to low-dimensional 
many-body systems only. 

The first step to achieve control over large-scale systems 
by controlling only a few particles was done in the context of 
quantum state transfer in spin chains. In these simple mod- 
els, the issues mentioned above were avoided by restricting 
the analysis to the subspace with a single excitation. Most 
of the proposed schemes for quantum state transfer are 
not actually based on the framework of control theory, but on 
smart tricks from various fields of physics, often using classi- 
cal intuition about the dynamics. Also, the theory of optimal 
control was recently applied to state transfer |7, 8]. 

A few of these proposals were then applied to more general 
tasks than state transfer; that is, how to use spin chains for en- 
tanglement purification [9], cloning transformations [10], or 
even for fully fledged quantum computation with little con- 
trol 1 6]. Thus, such schemes are no longer restricted to a 
small subspace of the full Hilbert space. These schemes for 
quantum computation use clever methods to design the con- 
trol pulses analytically, i.e., without relying on control theory. 
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but at the price of limitations on their appHcability. In par- 
ticular, they assume specific coupling parameters and design 
of the system. Nonetheless, such methods are intriguing from 
a theoretical perspective (e.g., in relation to complexity ques- 
tions) and give hope for the feasibility of quantum control of 
larger systems. 

The goal of this paper is to efficiently compute control 
pulses for a large system, using the full Hilbert space, and 
to show that the duration of the pulses scales efficiently (i.e., 
polynomially) with the system size. There are various trade- 
offs to consider For instance, if we allow more means for 
external control and/or higher ability in designing the system 
Hamiltonian, the problem will become more tractable theoret- 
ically (on paper), but more demanding and less relevant from 
an experimental perspective. 

Here, we present a solution to achieve feasible control both 
theoretically and experimentally. We will use a Hamiltonian 
that can be efficiently diagonalized for large systems through 
the Jordan-Wigner transformation. The control pulses are ap- 
plied only to the first two spins of a chain (see Fig. [U. The 
control consists of two parts: one where we will use the 
Jordan-Wigner to efficiently compute and control the infor- 
mation transfer through the chain (thus using it as a quantum 
data bus), and a second one where we will use some local 
gates acting on the chain end to implement two-qubit oper- 
ations. To be efficiently computable, these local gates need 
to be fast with respect to the natural dynamics of the chain. 
Combining the two actions allows us to implement any uni- 
tary operation described in the gate model. The interplay of 
the Jordan-Wigner picture for the information transfer and the 
gate model in the canonical basis will also require some tricks 
to avoid the accumulation of (uncomputable) phases. 

System and information transfer problem. — We consider 
a chain of N spin- 1/2 particles coupled by the Hamiltonian 
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c„ [(1 + 'y)XX + (1 - 'y)YY] 



where X, Y, Z are the Pauli matrices, the c„ are generic cou- 
pling constants and B„ represent a magnetic field. Varying 
the parameter 7 allows to encompass a wide range of exper- 
imentally relevant Hamiltonians 111 ll Il2l 11311 . including the 
transverse Ising model (7 = 1; for this case we require the 
fields Bn 7^ 0) and the XX model (7 = 0). We assume 
that the value of the first magnetic field, Bi , can be controlled 
externally. This control will be used to induce information 
transfer on the chain and realize swap gates 12Z] between ar- 
bitrary spins and the two 'control' spins 1, 2 at one chain end. 
Hence such swap gates are steered indirectly by only acting 
on the first qubit. The tunable parameter will also be used for 
applying unitaries on the control spins. 

In order to focus on the main idea we now present our 
method for 7 = and _B„ — for n > 1. The gen- 
eral case follows along the same lines, though being more 
technically involved. Our first task is to show that by only 
tuning Bi{t), we can perform swap gates between arbitrary 



pairs of qubits. First we rewrite the Hamiltonian, using 
the Jordan-Wigner transformation a„ = cr+ IlnKri ^m, into 
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^an}- The operators a„ obey 



the canonical commutation relations {a„,aT^} = Snm and 
{a„, a„i} = 0. The term we control by modulating Bi{t) is 
— Zi — 1 — 2a\ai. From the theory of quantum control 
^ we know that the reachable set of unitary time-evolution 
operators on the chain can be obtained from computing the 
dynamical Lie algebra generated by ihi and iH. It contains 
all possible commutators of these operators, of any order, and 
their real linear combinations. For example, it contains the 
anti-Hermitian operators i/ii2 =— [ihi,[ihi,iH]]/ {4ci) = 
i{a\a2 + o,\ai), ihi^ = [iH,ihi2] /c2 = 0.1O.3 — a\ai and 
ih23 = [ihi2,ihi3] = 1(0303 + 0302). We observe that 
taking the commutator with hi2 exchanges the index 1 of 
hi3 with 2. Taking the commutator with iH we find that 
ihi4 = [ihi3,iH] + ici/123 — ic2/ii2 = i{a\a4 + ajai) 
and ih24 = 0204 — a\a2 are also elements of the dynami- 
cal Lie algebra. Hence the effect of taking the commutator 
with H is raising the index of the hki ■ Generalizing this, we 
find that the algebra contains the elements ih^i, with k < I, 
ihki = alai - ajak for (fc - I) even, ihki = «(a|,a; + ajak) 
for {k — I) odd, and hk = Zk = 1 ~ 2aj,afe. From con- 
trol theory we thus know that the time evolution oper- 
ators exp{—3TTihki/2) (which will turn out to be very sim- 
ilar to swap gates) can be achieved through tuning Bi(t). 
The main point is that because both hi and H are free- 
Fermion Hamiltonians, the corresponding control functions 
can be computed efficiently in a 2A^-dimensional space (we 
will do so explicitly later). Ultimately, we need to transform 
the operators back to the canonical spin representation. Us- 
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for (fc — /) odd. The operator Lki = Y[k<j<i '^j ^ises from 
the non-local tail of the Jordan-Wigner transformation and act 
only on the state of the spins between k and I. They are con- 
trolled by the state of the qubits k,j being in the odd parity 
sector. 

In order to use the chain as a quantum data bus, our goal 
is to implement swap gates Ski = |00)fe;(00| + + 
10)fci(01| + |01)fci(10|,so the fact that we have achieved some 
modified operators with different phases on k,l instead, and 
also the controlled non-local phases Lki, could potentially be 
worrisome. It is worth pointing out that the XX interaction is 
not capable of generating swap gates of the (physical) qubits. 
We will use a method suggested in [ 14] that allows us to tackle 
these complications. That is, rather than using the physical 
qubits, we encode in logical qubits, consisting of two neigh- 
bouring physical qubits each. They are encoded in the odd 
parity subspace |01), |10). Swapping a logical qubit n to the 
control end of the chain then consists of two physical swaps 
cxp(— Tri/ii 2n-i/2) and cxp(— 7ri/i2 2ri/2). Since both phys- 
ical swaps give the same phases, the resulting operation is in- 
deed a full logical swap. Any single-qubit operation on the 
logical qubits can be implemented by bringing the target qubit 
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to the control end, performing the gate there, and bringing it 
back again. If we modulate the magnetic field Bi (t), we could 
equally decide to perform single logical qubit gates directly, 
without bringing them to the control end. This is possible be- 
cause exp(— i/i2n-i 2rii) in the physical picture translates to 
exp {—iXL_nt) in the logical picture, and because Z2n-i is in 
the algebra generated by Zi, which allows us to perform the 
operation exp {—iZ2n-it) = exp {—iZL^nt)- 

For quantum computation, we need to be able to perform 
at least one entangling two-qubit operation. We choose a 
controlled-Z operation, as this can be performed by operat- 
ing only on half of two logical qubits. That is, to perform 
a controlled-Z between logical qubit n and to, we bring the 
physical qubits 2n — 1 and 2m — 1 to the control end, per- 
form a controlled-Z between them, and bring them back. It 
is easy to check that again all unwanted phases cancel out. 
The controlled-Z could not be efficiently computed in the in- 
terplay with the many-body Hamiltonian H, because it cannot 
be generated by a quadratic Hamiltonian in the Jordan-Wigner 
picture. Therefore, this gate has to be implemented on a time- 
scale tg much faster than the natural evolution of the chain, 
i.e., tg <^ minj{l/cj}. We can soften this requirement by us- 
ing control theory to generate exp {—iZiX2t) by modulating 
/3i(i)Yi (this is a linear term in the Jordan-Wigner picture), 
and then using a fast Hadamard gate on the second site to ob- 
tain exp {—iZxZ-it), which, together with exp {—iZit) and 
exp {—iZ2t), gives the controlled-Z gate. This gives us quite 
a remarkable conclusion: except for a fast Hadamard on a sin- 
gle qubit, all other controls required for quantum computation 
can be computed efficiently in the Jordan-Wigner picture. 

Efficiency. — The crucial question left open in the above 
is: how long does it actually take to implement the gates? Un- 
fortunately, the theory of quantum control does not provide a 
general answer, though some interesting progress was recently 
reported 1 19]. In order to evaluate the efficiency, we have nu- 
merically simulated a range of chain lengths and studied the 
scaling of the logical swap operation time T with the (physi- 
cal) chain length N. For the control, we choose the local mag- 
netic field Bi (t) at the first site. We set the coupling strength 
constant, namely c„ = J V n. To provide evidence of a poly- 
nomial scaling, we set the simulation time = {N — 1)^, 
(all times are given in units of 1/ J and h — 1) and verify for 
each N that we can find a specific Bl{t) that performs the 
logical swap f23\. We quantify our success by calculating the 
error of the operation e = 1 — F, where F — {\tTU'^Ug\/N)'^ 
is the fidelity between the time evolution U and the goal uni- 
tary Ug. We find the function Bi{t) using techniques from 
optimal control theory Briefly, the procedure is as 

follows: (1) an initial guess is made for the function Bi{t); 
(2) we run the optimal control algorithm to generate a new 
Bi (t) which decreases the error of our operation; (3) steps 1 
and 2 are iterated until the final error reaches a defined thresh- 
old e. In practice, it suffices to choose a threshold which is 
of the same order of magnitude as the error introduced by the 
Hadamard gate. 

If the algorithm converges for each N and the correspond- 
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Figure 2: (color online) The optimized function (t) that produces 
the physical swap between n — 1 and n — 29 for the chain A'^ = 30 
in a time Tgo (in units of 1/J). The inset shows a zoom to confirm 
that the time-scale of the oscillations of BI (t) are slower than the 
dynamics of the chain. 



ing Ttv, giving the optimal pulse sequence -B*(t), then we 
can assert that the scaling of the operation time is at least 
as good as = {N — 1)-^, up to a given precision. We 
initialize the algorithm with a guess for the function Bi{t), 
e.g., Bi{t) = 1. This is, of course, a rather poor choice for 
performing the swap gate. However, the optimization algo- 
rithm exponentially improves the choice of Bi{t) reaching, 
after several thousand iterations, the desired precision. A typ- 
ical optimized function B^ {t) is given in Figure|2] It is worth 
pointing out that although the pulse looks like it is oscillating 
very rapidly, its oscillations are slower than the natural dy- 
namics of the chain, thus the control function does not need 
to be fast. Indeed, a Fourier analysis confirmed that only fre- 
quencies of up to ^ 0.02 J are required. Furthermore, such 
pulses are robust against small fluctuations 10]. The results of 
the optimizations for different chain lengths are shown in Fig- 
ure[3] the desired quadratic scaling has been clearly achieved. 
We stress here that the chosen scaling law may not neces- 
sarily describe the shortest time on which the physical swap 
gate can be performed. However, the dynamical Lie algebra of 
quasi-free fermions has a dimension of the order N'^, indicat- 
ing that such scaling might be optimal. Note that even though 
here we focused on swapping between the chain ends, which 
one expects to takes the longest time, we checked that general 
swaps between any pair of logical qubits can be achieved in 
(at least) the same time-scale (data not shown). 

A final remark on the robustness against imperfections of 
the results presented is needed, as perfect homogeneity or 
even fine tuning of individual couplings might be very hard to 
achieve in most systems and might prevent information trans- 
fer schemes to work [16]. In the previous analysis, for sim- 
plicity, we considered a chain with uniform couplings c„ — J 
for n > 1, but the results for arbitrary couplings are similar 
as long as localization effects can be neglected. We assume 
that the disordered Hamiltonian parameters are known, be- 
cause they can be estimated efficiently by controlling only the 
chain end BlTlllSll . We checked that the optimization with off- 
site disorder of 10%, uniformly distributed, leads to the same 
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Figure 3: (color online) The swap operation time Tjv versus the chain 
length A'^: the red continuous line is (TV — 1)^ while the black dots 
are the lengths for which we numerically verified this scaling with an 
error e < 10~^. 



results,for most realizations, until at least N — 40. 

Conclusion. — We have shown how to efficiently com- 
pute control pulses for large spin chains described by a vast 
class of Hamiltonians. The pulses are computed for an 
2Af-dimensional system but can be applied to the full 2^- 
dimensional system. Full quantum computation is possible by 
only controlling two spins at one end of the chain. The only 
price for this indirect control is that the quantum computation 
takes quadratically longer than for direct control. Given the 
large benefit of requiring so little control for a quantum com- 
puter, we believe that this scheme would be very useful for 
future implementations. As a further application, we remark 
that our proposal can also be applied to use the spin chain as 
a quantum memory, storing qubits by moving states from the 
controlled part to the rest of the chain and by applying con- 
trol pulses on qubits 1 , 2 to achieve the identity on the regis- 
ter (effectively switching off the chain Hamiltonian). For fu- 
ture studies, we would additionally like to probe the ultimate 
limit at which one may perform swaps using optimal control 
nil, Hill], and try to obtain simple (possibly analytic iflll ) 
pulses. 

As a possible application of our proposal, we note that in 
Josephson qubit implementations, much progress has been re- 
ported on the control of two-qubit gates fl 1], but that control- 
ling and reading out many qubits is very hard. In fact, it is 
expected to be difficult to construct fully addressable long ar- 
rays of Josephson junctions while it is conceivable to produce 
long chains of qubits with always-on interaction where only 
one or two qubits are fully controllable and readable 124)] . In 
those systems the decoherence time is (optimistically) ^ 1000 
times larger than the time-scale of the inter-qubit coupling 
ll20l,l2l]l . which would make our scheme applicable for up to 
^ 30 qubits to achieve a single swap gate (though, of course, 
much less to achieve full computation). Hence, the optimal 
control ideas presented here have the potential to address a se- 
rious limitation in such implementations and thus open a novel 
avenue towards quantum information processing in solid state 
devices. 
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